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Abstract 

o 

We study the rare semileptonic and radiative leptonic B s decays in the universal 
\ extra dimension model. In this scenario, with a single extra dimension, there exists 

only one new parameter beyond those of the standard model, which is the inverse of 
the compactification radius R. We find that with the additional contributions due 
to the KK modes the branching ratios of the rare B s decays are enhanced from their 



corresponding standard model values and the zero point of the forward backward asym- 
metries are shifted towards the left. 
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1 Introduction 



Although the standard model (SM) of electroweak interaction is very successful in explaining 
the observed data so far, but still it is believed that there must exist some new physics 
beyond the SM, whose true nature is not yet well-known. Therefore, intensive search for 
physics beyond the SM is now being performed in various areas of particle physics. In 
this context the B system can also be used as a complementary probe for new physics. 
Unfortunately, we have not been able to see any clear indication of physics beyond the SM 
in the currently running £?-factories (SLAC and KEK). Nevertheless, there appears to be 
some kind of deviation in some b — > s penguin induced transitions like the deviation in 
the measurement of sin2/3 in B d — > <pK s and also in some related processes, polarization 
anomaly in B — > <pK* and deviation of branching ratios from the SM expectation in some 
rare B decays, etc. [1]. But, the present situation is that it is too early to substantiate or 
rule out the existence of new physics in the b-sector. 

One of the important ways to look for new physics in the b-sector is the analysis of rare B 
decay modes, which are induced by the flavor changing neutral current (FCNC) transitions. 
The FCNC transitions generally arise at the loop level in the SM, and thus provide an 
excellent testing ground for new physics. Therefore, it is very important to study the FCNC 
processes, both theoretically and experimentally, as these decays can provide a sensitive test 
for the investigation of the gauge structure of the SM at the loop level. Concerning the 
semileptonic B decays, B — > X s l + l~ (X s = K,K*, I = e,fj,,r) are a class of decays having 
both theoretical and experimental importance. At the quark level, these decays proceed 
through the FCNC transition b — > sl + l~, which occur only through loops in the SM, and 
therefore, these decays constitute a quite suitable tool of looking for new physics. Moreover, 
the dileptons present in these processes allow us to formulate many observables which can 
serve as a testing ground to decipher the presence of new physics. In this context, the study 
of the rare semileptonic B s decay modes, B s — > which are induced by the same quark 

level transition, i.e., b — > sl + l~ , might be worth exploring. These decay modes may provide 
us additional information towards the quest for the existence of new physics beyond the SM 
and therefore deserve serious attention, both theoretically and experimentally. Since at the 
quark level they are induced by the same mechanism, so we can also independently test our 
understanding of the quark-hadron dynamics and also study some CP violation parameters 
with the help of rare B s decays, apart from corroborating the finding of the 5-meson sector. 
These decay modes are studied in the standard model [2] and in the two Higgs doublet 
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model [3]. Recently, the DO collaboration [4] has reported a more stringent upper limit on 
the branching ratio of B® — > mode as Br(£?° — > < 4.1 x 1CT 6 . 

Another sensitive process to look for new physics is the radiative dileptonic decay mode 
— > l + l^ , ~f. In contrast to the pure leptonic decays B° s — > l + l~ , these modes are free from 
helicity suppression due to the emission of a photon in addition to the lepton pair. Thus, the 
branching ratios of these processes are larger than those for the pure leptonic modes despite 
of an additional a suppression. These modes are also studied in the SM [5] and some beyond 
the SM scenarios [6]. 

New physics effects manifest themselves in these rare B s decays in two different ways, 
either through new contribution to the Wilson coefficients or through the new structure in 
the effective Hamiltonian, which are absent in the SM. There are many variants of possible 
extensions to the SM exist in the literature but the models with extra dimensions have 
received considerable attention in recent years. An elegant beyond the SM scenario with 
extra dimensions is being considered to solve the classical problems of particle physics, mainly 
the hierarchy problem, gauge coupling unification, neutrino mass generation, fermion mass 
hierarchies etc. In fact, there are several models that exist in the literature which take 
into account the effect of extra dimensions and they differ from one another depending on 
the number of extra dimensions, the geometry of space-time, the compactification manifold, 
which particles can go into the extra-dimensions and which cannot etc. In the case of the 
scenario with Universal Extra Dimensions (UED), proposed by Appelquist-Cheng-Dobrescu 
(ACD) [7], all the fields are allowed to propagate in all available dimensions. In its simplest 
version the single extra dimension is taken to be the fifth dimension, that is, £5 = y. This 
is in fact compactified on the orbifold S 1 /Z 2 , i.e., on a circle of radius R and runs from 
to 27iR, with y = and y = ttR are the fixed points of the orbifold. Hence a field F(x,y), 
where x denotes the usual 3 + 1 dimension, would be a periodic function of y and can be 
represented as 

00 

F(x,y)= £ F n (x)e in -y. (1) 

n=— oo 

Under the parity transformation P5 : y —y, fields which exist in the SM are even and their 
zero modes in the KK expansion are interpreted as ordinary SM fields. On the other hand, 
fields absent in the SM are odd under P 5 , so they do not have zero modes. Compactification of 
the extra dimension leads to the appearance of Kaluza-Klein (KK) partners of the SM fields 
as well as KK modes without corresponding SM partners. Thus, in the particle spectrum of 
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ACD model, in addition to the ordinary particles of the SM denoted as zero mode (n — 0), 
there are infinite towers of Kaluza-Klein modes (n > 1). There is one such tower of each SM 
boson, two for each SM fermion, and for the physical scalars a ^, no zero mode exists. The 
masses of the bosonic KK modes are given as m 2 n ^ = ml + n 2 /R 2 , with m being the mass of 
the zero mode. The important features of the ACD model are: (i) there is only one additional 
free parameter with respect to the SM, i.e., l/R, the inverse of the compactification radius 
(ii) the KK parity is conserved, which implies that the KK modes do not contribute at tree 
level for the low energy interaction processes and (iii) the lightest KK particle must be stable. 

The implications of physics with UED are being examined with the data from accelerator 
experiments, for example, from Tevatron experiments the bound on the inverse compacti- 
fication radius is found to be about l/R > 300 GeV. Analysis of the anomalous magnetic 
moment and Z — ► bb vertex [8] also lead to the bound l/R > 300 GeV. Possible manifes- 
tation of ACD model in the Kl — K$ mass difference, B° — B° mixing, rare decays of K 
and B mesons and the kaon CP violation parameter e'/e are comprehensively investigated 
in [9]. Exclusive B — > K*£ + £~, B — > K*vv and B — > K*^y decays [10] and rare semileptonic 
A5 decays [11] are also studied in the framework of the UED scenario. 

In this work, we would like to study the implications of the scenario with universal extra 
dimensions in the semileptonic rare decay modes B° s — > <fi , and the radiative dileptonic 
decays B° s — > As discussed above, since these decays are very sensitive to new physics, 

these studies may provide an indirect way to constrain the new physics parameter R. The 
paper is organized as follows. In section II, we present the effective Hamiltonian for these 
decays in the ACD model. The decay rates and forward backward rate asymmetries for the 
mode B® —> <fil + l~ and B® — > Z + Z~7 are considered in section III. Section IV contains our 
conclusion. 



2 Effective Hamiltonian 



The effective Hamiltonian describing the b — > sl + l transition in the SM can be written as 
[12] 



G F a 



C^is^Lb^Yl) 



+ C 10 (s^Lb)(h^ 5 l) ~ 2C e 7 ff m b (sta^Rb)(iri) 



(2) 
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where q is the momentum transferred to the lepton pair, given as q = p- + p+, with p_ 
and p + are the momenta of the leptons l~ and / + respectively. L,R = (1 ± 7s)/2 and Cj's 
are the Wilson coefficients evaluated at the 6 quark mass scale. The coefficient Cg^ has a 
perturbative part and a resonance part which comes from the long distance effects due to 
the conversion of the real cc into the lepton pair l + l~ . Hence, can be written as 

C e 9 ff = C 9 + Y(s) + Cl es , (3) 

where s = q 2 and the function Y(s) denotes the perturbative part coming from one loop 
matrix elements of the four quark operators and is given in Ref. [13]. The long distance 
resonance effect is given as [14] 

C r 9 es = ^C 1 + C 2 + 3C 3 + C, + 3C 5 + C 6 ) £ K ™v i W-+l + l-) 
a ViJjfw rn 2 Vi -s-im Vi Y Vi 

where the phenomenological parameter k is taken to be 2.3, so as to reproduce the correct 
branching ratio B(B -> J/^K* -> K*l+l~) = B(B -> J/ipK*)B(J/ip -> l+l~). 

Now we will analyze the KK contributions to the above decay modes. Since there is 
no tree level contribution of KK modes in low energy processes, the new contributions, 
therefore, only come from the loop diagrams with internal KK modes. It should be noted 
here that there does not appear any new operator in the ACD model, and therefore, new 
effects are implemented by modifying the Wilson coefficients existing in the SM, if we neglect 
the contributions of the scalar fields, which are indeed very small. Thus, the modified Wilson 
coefficients are given as [9] 

11 a s (M w ) 11 a s (M w ) 

C 3 (M W ) = -3C 4 (M w )=C 5 (M w ) = -3C 6 (M w ) = --fE(x t ,l/R), 

Z47T 

C 7 {M W ) = ~&(xt,l/R), C s {M w ) = - l -E'{x t MR). (5) 
The functions E(x t ,l/R), E'(x t , 1/R) and D(x t ,l/R) are given by 

oo 

F(x t , 1/R) = F (x t ) + F n(*t, x n ), F — E, E', D' , (6) 

71=1 

where x t = mf/M^, x n = m^/M^ and m n = n/R. The loop functions F (x t ) and F n (x t , x n ) 
are given as [9] 

, , 2, x?(15- l§xt + ±x 2 ) , x t (18 - llx t - x 2 ) 

E (x t ) = --lnx t + ^— t—^L\nx t + ^— , 

3 6(1 -x t ) A 12(1 -x t ) s 
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(8x* + 5x?-7x t ) x 2 t (2-3x t ) 

Do{Xt) ~ 12(1 - Xt f + 2(1 - Xt r " 

km = - Xt{ i: hx ::, 2) AtAtv^^ 



4(1 -x t f 2(1 -x t y 
x t [35 + 8a; t - 19x 2 + 6^(10 - 9x t + 3x 2 ) + 3x n (53 - 58x t + 21a; 2 )] 



36(x t - l) 3 



1 '1 + x„)(-2 + 3x n + x 2 n ) In ■ 



D|,(a; t ,x n ) = 



2 V " ' - ,v " ' ™" ' ~ ny ~"l + a; n 
(1 + x n )(-6 + lggf - 9x 2 + x 2 (3 + x t ) + x n (9 - 4x f + 3a: 2 )) x„ + x f 

6(a^-l) 4 11 l + x n ' 

x t [-37 + 44a; t + 17x 2 + 6x 2 (10 - 9a;t + 3x 2 ) - 3x„(21 - 54a; t + 17x 2 )] 



36(:r t - If 



Xn 

(2 — 7x n + 3x 2 ) ^ x 



1 + av, 



(-2 + a: w + 3xf)[x f + 3a: 2 + x 2 (3 + x t ) - x n (l + (-10 + a; t )x f )] x TO + x t 

6(a; t -l) 4 11 l + x n ' 

x t [-17 - 8a; t + x\ - 3x n (21 - 6x t + a; 2 ) - 6a; 2 (10 - 9x t + 3a; 2 )] 
" I2(x t - l) 3 

- -x n (l + x„)(-l + 3a; n ) In — 
2 1 + x„ 

(1 + x TO )[a: f + 3a: 2 + a: 2 (3 + x t ) - x n (\ + (-10 + a; f )a: t )] a:„ + rr f 

2(x t - l) 4 11 1 + x n ' U 

The functions with and without x n correspond to the KK excitation and SM contributions, 
respectively. The summations are carried out using the prescription presented in [9]. 

After obtaining the Wilson coefficients at the scale M W) we have to run these coefficients 
Ci(M w ) down to the scale /j, ~ m b . The Wilson coefficients Cj=i ; ... 6 (// ~ m b ) at low energy 
can be obtained from the corresponding C i= i^... & (M w ) ones by using the Renormalization 
Group (RG) equation, as discussed in Ref. [12], as 

C(fi) = U 5 (fi,M w )C(M w ) , (8) 

where C is the 6x1 column vector of the Wilson coefficients and Us(yU, My/) is the five-flavor 
6x6 evolution matrix. In the next-to- leading order (NLO), Us(/x, Mw) is given by 

U 5 („, M W ) = (l + S^j) Uf („, M W ) 1 1 - ^-j) , (9) 
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where uf\n,M w ) is the leading order (LO) evolution matrix and J denotes the NLO 
corrections to the evolution. The explicit forms of U 5 (//, M w ) and J are given in Ref. [12]. 
For the coefficient C-j we use the RG running as 



WJ 



7/< 



(10) 



i=l 



where rj = ^j^y, and the coefficients ai and hi are as given in Ref. [9]. 
The Wilson coefficient Cg in the ACD model and in the NDR scheme is 

,NDR , Y(x t ,l/R) 



+ 



sin 2 6 



4Z(x t , l/R) + P E E(x t , l/R) , 



(11) 



w 



where P^ DR = 2.60 ± 0.25 [9] and the last term is numerically negligible. The functions 
Y(x t ,l/R) and Z(x t ,l/R) are defined as: 



Y(x t ,l/R) = Y (x t ) + J2C n (x t ,x n ) , 



n=l 

00 



Z(x t ,l/R) = Z (x t ) + C n (x t ,x n ) , 

n=l 



(12) 



with 



Yo(x t ) = I 



x t -4 



+ 



3^ 



^ - 1 (x t - 1) ; 



In 



18^ - 1634 + 2594 - 108x t 



+ 



UA(x t - l) 3 
- , - 38x 3 t - 15x 2 t 



32x\ - 38x 3 t - 15x 2 + lSx t 1 



72(x t - iy 



\nx t 



C n (xti x nj 



x t 



8(x t - If 



x\ — 8x t + 7 + (3 + 3x t + 7x n - x t x n ) In 



1 + x n 



(13) 
(14) 



Finally, the Wilson coefficient C w , which is scale independent, is given by 

Y(x t ,l/R) 



C 



10 



sin 2 9 W 



(15) 
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3 Decay rates and forward backward asymmetries 



3.1 B s — > (j) l + l process 



The semileptonic rare B decays are very well studied in the literature, both in the SM 
and in various extensions of it. Therefore, we shall simply present here the expressions 
for the dilepton mass spectrum and the forward backward asymmetry parameters. For the 
B s — > (f)l + l~ process, the hadronic matrix elements are given as 

(<Kk,e)\(V - A),\B S (P)) = e^e^P^^p- - te;(m B + mJA^q 2 ) 

m B + 



+ i(P + k) fl (e*q) 



Mq2) +^%{e* q ) 2 -^\M<f)-M ( ?) 



((f>(k,e)\sa lll/ q u (l + ~f 5 )b\B(P)) = ie^e** ■P a k f3 2T 1 (q 2 ) + 



* / 2 2 ^ 



-(e*q)(P + k), T 2 (q 2 ) + (e*q) 



m 2 B -m 2 



(P + k),^ 2 ) , 



(16) 



where V and A denote the vector and axial vector currents, A , A±, A 2 , A s , V, T±,T 2 and T 3 
are the relevant form factors and m B denotes the mass of B s meson. Thus, with eqs. (2) 
and (16) the transition amplitude for B s — > (j)l + l~ is given as 



M(B S 



= ^^V u yz{wi[2A £l ^e"k a pf' + iBe m lt 

- iC(P + k)^(e*q) - tD(e*q)q,} + lYl 5 l\2Ee^e*"k a P^ 



+ iFe;-iG(e*q)(P + k) lx -iH(e*q)q tl 
where the parameters A,B,---H are given as [3] 



(17) 



A 
B 

C 
D 
E 



C e 9 ff ^^- + 4^C^T 1 (q 2 ), 



lm B 



+ m<t> ) (c^A^q 2 ) + 4^(m B - m^C e 7 ff T, 



= a 



l"^p- + 4^C?" (T 2 (q 2 ) + ^— 
m B + 171^ q 2 \ rn%- ra£ 

2C^(A 3 (q 2 ) - A {q 2 )) - ACf^T 3 (q 2 ) , 

Cio V{q2) , F = C w (m B + m (f) )A 1 (q 2 ) , 
m B + 



2{q 



Ts(q 2 ) 



G 



C 



10" 



m B + 

The differential decay rate is given as 

dT GW 



, H = 2C w '-^(A 3 (q 2 )-A (q' 



(18) 



ds 2 12 ir 5 rriB 



|^ t ;| 2 A 1/2 (l,r^) vi A , 



(19) 



where s = q 2 /m 2 B , = m 2 ^/m 2 B , v\ = \Jl — Amf/s, A = A(l,r^,s), is the triangle function 
and 



A = -\m B s((3-vf)\A\ 2 + 2vf\E\ 2 ) + — Am 
3 v ' r> 



-\m 2 B (3 - v 2 )\C\ 



+ m 2 B s 2 (l - vf)\H\ 2 + - [(3 - ^ 2 )(r + s - 1) - 3s(l - v 2 )] Re[FG* 
- 2s(l - vf )Re[FH*} + 2m|s(l - t>)(1 - v 2 )i?e[Gi/*] 



+ 3(3-u, 2 )(r + s-l)ite[SC* 



+ 



3t\/ 



-m. 



(A + 12r>s)(3 - t; 2 )|5| 2 + Am* [A(3 - v 2 ) 



- 3s(s - 2r> - 2)(1 - u, 2 )] |Gf + (A(3 - vf) + 2Arfsv 2 ) \F\ 



(20) 



Another observable is the lepton forward backward asymmetry (A FB ), which is also a very 
powerful tool for looking new physics. The position of the zero value of Ap B is very sensitive 
to the presence of new physics. The normalized forward-backward asymmetry is defined as 



'-FB 



» 



i d 2 T 
o dsd cos 9 



dcos9 



° d 2 T 



dsd cos 9 



dcos9 



I 



d cos 9 ■ J — • - ■ —(la >s 9 



dsd cos 



dsd cos 9 



GW 



* |2 



2 12 vr 5 m s 



8ml; A w 2 s (Re[BE*\ + i?e[AF* 



(21) 



where 9 is the angle between the directions of l + and B s in the rest frame of the lepton pair. 

For numerical evaluation we use the form factors calculated in the LCSR [15] approach, 
where the q 2 dependence of various form factors are given by simple fits as 

T2 



2 ' 



1 - q 2 /m) it 

n 

1 - q 2 /m 2 R 1 - q 2 /mj it ' 



r 2 



ri 



/(g ) 1 - q 2 /m 2 



+ 



r 2 



fit 



(for A u T 2 ) 
(for V, Aq, 7i) 
(for A 2 , T 3 ) • 



(22) 
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The values of the parameters r±, r 2 , rriR and mf it are taken from [15]. The uncertainties 
occur in these fitted results are also given in [15]. The form factors A 3 and T 3 are given as 

a i i\ m B + m v 2 m B -m$ 2 

Mq) = ~^r Mq) — 2^ Mqh 

T 3 (g 2 ) = ^^(T3(g 2 )-T 2 (g 2 )). (23) 

The particle masses, lifetime of B s meson and the weak mixing angle as sin 2 9w are taken 
from [16]. The quark masses (in GeV) used are m b =4.6, m c =1.5 and the CKM matrix 
elements as V tb V t * s = 0.04 . 

In order to see the effect due to the uncertainties in the form factors, we plot the dif- 
ferential branching ratio (19) and the forward backward asymmetry (21) for B s — > <fi 
in Figure-1. It can be seen from the figure that due to the uncertainties in the form fac- 
tors these distributions deviate slightly from their corresponding central values in the low s 
region. However, in the large s region these deviations are highly suppressed and the zero 
position in the forward backward asymmetry is insensitive to these uncertainties. Therefore, 
we will not consider the effect of these uncertainties for the differential decay distributions 
and the forward backward asymmetries but we will incorporate their effects in the total 
decay rates. 




5 10 15 20 5 10 15 20 

s s 

Figure 1: The variation of the differential branching ratio dBi/ds (in units of 10~ 7 ) and the 
forward backward asymmetry in the SM with s (in GeV 2 ) for the process B s — > (f)fi + fi~ . 
The solid line denotes the result using the central values of the form factors and the dashed 
lines represent the effects due to the uncertainties in these parameters. 

Next, to see the effect of the UED, we plot the differential decay rates (19) and the 
forward backward asymmetries (21) for B s — > l + l~ against s, as depicted in Figures 2 and 
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3. It can be seen from these figures that there is considerable enhancement in the decay rates 
due to the KK contributions for 1/R — 200 GeV. The zero position of the forward backward 
asymmetry A FB shifts towards the left due to the NP effect. This shifting is also more 
prominent for 1/R = 200 GeV. Therefore, its experimental determination would constrain 
the parameter 1/R. 




s s 

Figure 2: The variation of the differential branching ratio dBr/ds (in units of 10~ 7 ) and the 
forward backward asymmetry with s (in GeV 2 ) for the process B s — > <j)n + fi~ . The solid line 
denotes the SM result, the dashed line represents the contribution from UED model with 
1/R = 200 GeV and the long-dashed line for l/i?=500 GeV. 




14 16 18 20 14 16 18 20 

s s 



Figure 3: Same as Figure-2 for the process B s — > 0t + t 

We now proceed to calculate the total decay rate for B s — > l + l~ . It should be noted 
that the long distance contributions arise from the real cc resonances with the dominant 
contributions coming from the low lying resonances J/ip and ip'. In order to minimize the 
hadronic uncertainties it is necessary to eliminate the backgrounds coming from the resonance 
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Figure 4: The variation of the total branching ratio in units of 1CT 7 with 1/R. The solid 
line represents the central value while the dashed lines are due to the uncertainties in the 
form factors. Figure-(a) represents the B s — > (f)fi + fi~ and figure (b) B s — > 0t + t _ process. 
The horizontal dot-dashed lines in both the figures represent the central values in the SM 
and the dashed lines are the corresponding uncertainties. 

regions. This can be done by using the following cuts for B s — > process as [17] 

2.9 GeV < m M+M - < 3.3 GeV , and 3.6 GeV < m M+M - < 3.8 GeV , 

and for B s — > 0t + t~ process m^> — 0.08 < m T + T - < m^i + 0.08 . 

Using these veto windows we obtain the branching ratios for the semileptonic rare B s 
decays which are given in Table-1. It is seen from the table that the branching ratios obtained 
in the ACD model are enhanced from the corresponding SM values for lower 1/R value. The 
variation of the total decay rates with 1/R are depicted in Figure-4. Thus, observation of 
these decay modes can be used to constrain the compactification radius R. As seen from the 
figure for large 1/R the standard model results are recovered. 

3.2 B° s -> Z+Z"7 process 

Now let us consider the radiative dileptonic decay modes B s — > / + /~7, which are also very 
sensitive to the existence of new physics beyond the SM. Due to the presence of the photon 
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Table 1: The branching ratios for various rare semileptonic (in units of 1CT 7 ) and radiative 
leptonic (in units of 1CT 9 ) B s decays, where Br SM represents the SM branching ratio and 
Br|i/R =2 oo is the branching ratio with KK contributions for 1/R = 200 GeV and Br^/R^oo 
for 1/R = 500 GeV. The errors are due to the uncertainties in the form factors. 



Decay process 


Br SM 


Br i/r=2oo 


Br i/r=5oo 


Expt.[4] 


B s — >■ 


15.29 ±0.95 


20.82 ± 1.06 


16.01± 0.99 


< 41 


B s -> (f)T+T- 


0.87±0.02 


1.53 ± 0.03 


0.98±0.025 




B s -> /i + /i"7 


3.05±0.50 


4.45± 0.69 


3.27±0.53 




B s — > r + r _ 7 


13.93±0.09 


24.93±0.14 


15.77 ± 0.10 





in the final state, these decay modes are free from helicity suppression, but they are further 
suppressed by a factor of a. However, in spite of this a suppression, the radiative leptonic 
decays B s — > / + /~7, / = (//, r) have comparable decay rates to that of purely leptonic ones. 

The matrix element for the decay B s — > Z + Z - 7 can be obtained from that of the B s — > l + l~ 
one by attaching the photon line to any of the charged external fermion lines. In order to 
calculate the amplitude, when the photon is radiated from the initial fermions (structure 
dependent (SD) part), we need to evaluate the matrix elements of the quark currents present 
in (2) between the emitted photon and the initial B s meson. These matrix elements can be 
obtained by considering the transition of a B s meson to a virtual photon with momentum 
k. In this case the form factors depend on two variables, i.e., k 2 (the photon virtuality) and 
the square of momentum transfer q 2 = (p B — k) 2 . By imposing gauge invariance, one can 
obtain several relations among the form factors at k 2 = 0. These relations can be used to 
reduce the number of independent form factors for the transition of the B s meson to a real 
photon. Thus, the matrix elements for B s — > 7 transition, induced by vector, axial-vector, 
tensor and pseudo-tensor currents can be parametrized as [18] 

F A 



(7(k,e)\s^ 5 b\B s (p B )) 

( 7 (M)|s7^M> 

(l(k,e)\sa IJ , l/ q l/ -f 5 b\B s (p B )) 
{l{k,e)\sa llu q v b\B s {p B )) 



^(Pb ■ k) - (e* ■p B )k fl 

k^, 
m B 

^.{Pb ■ k) - (e* ■p B )k^ 
k?F TV , 



ze 



m B 



(24) 
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where e and k are the polarization vector and the four-momentum of photon, p B is the 
momentum of initial B s meson and Fj's are the various form factors. 

Thus, the matrix element describing the SD part takes the form 

M SD = ^^Vuy^^e^plk^Aj^l + A^l^) 

+ t(e;(k-p B )-(e*-p B )k tl )(B 1 lYl + B 2 lYl 5 l)} , (25) 

where 

A\ = 2Cj—F TV + C 9 — — , A 2 = Cio — — , 

q 2 m B m B 

Bi = —ICj—^-Fta — C%—^- , B 2 = — Ciq—^- . (26) 

q l m B m B 

The form factors F v and F A have been calculated within the dispersion approach [19]. The 
q 2 dependence of the form factors are given as [18] 

HE,) = f>£f§r , (27) 
where E 1 is the photon energy, which is related to the momentum transfer q 2 as 

^M'-i)- (28) 

The values of the parameters (3 and A are given in Table-2. The same ansatz (27) has also 
been assumed for the form factors F TA and F TV . The decay constant of the B s meson is 
not yet well known because the pure leptonic decays of B s meson (i.e., B s — > l + l~) from 
which it could be extracted are highly suppressed in the SM as they occur at one one-loop 
level. Using QCD Sum rule approach, its value is found to be f Ba = 236 ± 30 MeV [20] 
and f Bs = 244 ± 21 MeV [21]. The Lattice QCD result gives f Ba = 242 ± 9 ± 34 MeV [22]. 
Therefore, in this analysis we use the value f Bs = 240 MeV. 

When the photon is radiated from the outgoing lepton pairs, the internal bremsstrahlung 
(IB) part, the matrix element is given as 



a 3 / 2 G F 

Mib = — -j=- V tb V* f Bs m t do 



n— h -— h w 

. \p+ ■ k p- ■ k J 



(29) 



Thus, the total matrix element for the B s — > l + l 7 process is given as 

M = M SD + M IB . (30) 
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Table 2: The parameters for B — > 7 form factors. 



Parameter 




Ftv 


F A 


Ft a 


^(GeV- 1 ) 
A(GeV) 


0.28 
0.04 


0.30 
0.04 


0.26 
0.30 


0.33 
0.30 



The differential decay width of the B — > l + l 7 process, in the rest frame of B s meson is 
given as 

dT G%a 3 2 3 

™s Ai , (31) 



ds 2 10 7T 4 



where 



Ai = ^m|(l-s) 2 ^((s + 2r i )(|A 1 | 2 + | J B 1 | 2 ) + (s-4r i )(|A 2 | 2 +|5 2 | 2 ) 

- 64 AT^^o((4n-S 2 -l)lnl±^ + 2S^) 

m^ s 1 — s v 1 — vi ' 

- 32n(l-s) 2 f Bs Re(C w Al), (32) 



with s = q 2 , s = s/m 2 B , T\ = mf/m 2 B , v\ = \Jl — Amf/q 2 . The physical region of s is 
4m 2 < s < m 2 B . 

The forward backward asymmetry is given as 

1 



Afb ~ AT 



2m 2 B s{\ - sfvf Re(A\B 2 + B{A 2 ) 



+ 32 f Bs n(l - s) 2 \n (^±)Re(c 10 B* 2 ) 



V s 



(33) 



Again to visualize the effects due to the uncertainties in the form factors (assuming it to 
at the level of 10%) we first plot the differential decay distribution (31), and the forward 
backward asymmetry (33) for B s — > yU + /i _ 7 in Figure-5. From the figure it can be seen that 
these uncertainties can affect only the differential decay rate but not the forward backward 
asymmetry. Furthermore, it is found that these effects are significant only for B s — > yU + yU~7 
and not for B s — > t + t~7. Next, we would like to see the effect of UED in these differential 
distributions. For this purpose, we consider only the central values of the form factors, 
/b s =0.24 GeV, a=l/128 and plot plot the dilepton mass spectrum (31), and the forward 
backward asymmetries (33) for B s — > Z + /~7 decays which are shown in Figures-6 and 7. From 
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these figures, we see that the branching ratio for B s — > Z + /~7 enhanced significantly from their 
corresponding SM values. However, the forward backward asymmetry is reduced slightly 
from the corresponding SM value for the B s — > process and there is a backward 

shifting of the zero position. For B s — > t + t~7 process there is no significant change in the 
asymmetry distribution. 




5 10 15 20 25 30 5 10 15 20 25 30 

s s 

Figure 5: The differential branching ratio (31) (in units of 10~ 9 ) and the forward backward 
asymmetry (A FB ) for the process B s — > /i + /i"7 in the standard model where the solid line 
represents the central value and the dashed lines represent the uncertainties due to the form 
factors. 

To obtain the branching ratios it is necessary to eliminate the background due to the 
resonances J/t/j(ip') with J/t/j(ip') — > l + l~ . We use the following veto windows to eliminate 
these backgrounds 

B s ->• /i + /i~7 : m J/4> ~ 0-02 < m M+M - < m J/v , + 0.02; 

: m^i — 0.02 < m M + M - < m^> + 0.02 
B s — > r + r _ 7 : m^i — 0.02 < m T + T - < m,^ + 0.02 . 

Furthermore, it should be noted that the |-M/,b| 2 has infrared singularity due to the emission 
of soft photon. Therefore, to obtain the branching ratio, we impose a cut on the photon 
energy, which will correspond to the experimental cut imposed on the minimum energy 
for the detectable photon. Requiring the photon energy to be larger than 25 MeV, i.e., 
E,y>5 m Bs /2 ) which corresponds to s < m Bs (l — 5), and therefore, we set the cut 5 > 0.01. 

Thus, with the above defined veto windows and the infrared cutoff parameter, we obtain 
the branching ratios as shown in Table-1 which are enhanced from their SM values. It should 
be mentioned that the B° s t + t~^ could be observable in the Run II of Tevatron. 
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Figure 6: The differential branching ratio (31) (in units of 10~ 9 ) and the forward backward 
asymmetry (A FB ) for the process B s — > yU + /i~7 in the standard model and in the ACD 
model. The solid line denotes the SM result, the dashed line represents the contribution 
from UED model with 1 /R = 200 GeV. 

4 Conclusion 



It is now widely believed that the scenario with extra dimensions is a strong contender to 
reveal physics beyond the SM and has received considerable attention in the literature. In 
view of this anticipation it is also worthwhile to study its implications in the b-sector. In this 
paper, therefore, we have studied the rare semileptonic decay mode B s — > <fi l + l~ and the 
radiative leptonic B s — > / + /~7 in the model with a single extra dimension. The branching 
ratios for various decay modes are found to be larger than their corresponding SM values 
but they lie within the present experimental upper limit. Furthermore, the zero point of the 
forward backward asymmetries for the decays under consideration are shifted to the left and 
the change is found to be sensitive to the inverse compactification radius. In future, with 
more intensive data, the UED scenario will be subjected to more stringent tests and in turn 
will enrich us with a better understanding of the flavor sector. 
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